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Abstract 


Prognostic cloud schemes are increasingly used in weather and climate models in 
order to better treat cloud-radiation processes. Simplifications are often made in such 
schemes for computational efficiency, like the scheme being used in the National Centers 
for Environment Prediction models that excludes some microphysical processes and 
precipitation-radiation interaction. In this study, sensitivity tests with a 2D cloud 
resolving model are carried out to examine effects of the excluded microphysical 
processes and precipitation-radiation interaction on tropical thermodynamics and cloud 
properties. The model is integrated for 10 days with the imposed vertical velocity denved 
from the Tropical Ocean Global Atmosphere Coupled Ocean— Atmosphere Response 
Experiment. The experiment excluding the depositional growth of snow from cloud ice 
shows anomalous growth of cloud ice and more than 20% increase of fractional cloud 
cover, indicating that the lack of the depositional snow growth causes unrealistically 
large mixing ratio of cloud ice. The experiment excluding the precipitation-radiation 
interaction displays a significant cooling and drying bias. The analysis of heat and 
moisture budgets shows that the simulation without the interaction produces more stable 
upper troposphere and more unstable mid and lower troposphere than does the simulation 
with the interaction. Thus, the suppressed growth of ice clouds in upper troposphere and 
stronger radiative cooling in mid and lower troposphere are responsible for the cooling 
bias, and less evaporation of rain associated with the large-scale subsidence induces the 
drying in mid and lower troposphere. 
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Introduction 


Clouds microphysical processes and cloud-radiation interaction play an important 
role in tropical weather and climate. Convective processes, which generate precipitation, 
cannot be resolved in general circulation models (GCMs). To meet the needs of the 
models, various cumulus parameterization schemes have been designed to estimate the 
amount of the precipitation and the associated heating and moistening effects using the 
large-scale variables in the GCMs. The most commonly used cumulus parameterization 
schemes in the GCMs have led to many successes in the simulation of atmospheric 
circulations. However, most of the schemes calculate precipitation diagnostically and do 
not allow an explicit representation of clouds and the cloud-radiation interaction in the 
GCMs. 

Sundqvist (1978) and Sundqvist et ai. (1989) first developed a prognostic cloud 
scheme for cloud water. The scheme allows an explicit representation of non- 
precipitating water clouds and the exchange of cloud water among the columns. They 
showed that the large-scale and mesoscale models with the scheme yield reasonable time 
evolution of cloud water and precipitation amounts. The prognostic cloud schemes have 
been applied to the GCM and regional models that led to improvement of numerical 
weather prediction in global and regional scales (e.g., Golding, 1990; Smith, 1990; 
Pudykiewicz et al., 1992; Tiedtke, 1993; Mannoji, 1995; Willoughby et al., 1984; Lord et 
al., 1984; Liu et al., 1997; Zhao and Carr, 1997). 

Among the many studies discussed above, Zhao and Carr (1997), hereafter 
referred to as ZC97, developed a prognostic cloud scheme with the diagnostic 
precipitation for operational numerical weather prediction models. The scheme includes 



prognostic equations for non-precipitating water and ice clouds and sophisticated 
microphysical processes associated with the growth of clouds and generation of 
precipitation. They showed that the inclusion of ice clouds in the models is important in 
transferring water vapor to precipitation and in the enhancement of latent heat release, 
which is consistent with the results by Willoughby et al. (1984), Lord et al. (1984), and 
Liu et al. (1997) that the inclusion of ice microphysics parameterization schemes in the 
models led to more realistic simulation of cloud structures of model hurricanes. The 
scheme has been used in the operational Global Data Assimilation Systerm (GDAS) at 
National Oceanic and Atmospheric Administration (NOAA)/NationaI Centers for 
Environmental Prediction (NCEP), USA. Since it is designed as a computational efficient 
scheme for operational models, ZC97 excludes some microphysical processes and 
calculates precipitation diagnostically without the physical presence of corresponding 
hydrometeors, which is equivalent to exclude precipitation-radiation interaction. 

In this study, tropical convective responses to microphysical and radiative 
processes are investigated with a 2D cloud resolving model to examine the possible 
effects of the microphysical processes and precipitation-radiation interaction on tropical 
thermodynamics and cloud properties excluded by ZC97. Concerning the fact that the 
ZC97 scheme and the bulk cloud microphysics scheme are not the same, and that the two 
schemes are applied to different models, the current study is more a sensitivity test rather 
than an evaluation of the effect of certain neglected microphysical processes. In section 2, 
the model forcing, and experiment designs will be discussed. In section 3, the model 
sensitivity tests and the comparison between the satellite-retrieved cloud products and 
GDAS data are carried out to show that the exclusion of the microphysical processes by 
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ZC97 in the model physics may simulate unrealistically large mixin g ratio of cloud ice. ■ 
The cause of the unrealistic simulation of cloud ice is further examined analyzing the 
budget of cloud ice in the cloud resolving simulations. In section 4, effects of 
precipitation-radiation interaction on the tropical thermodynamics are examined using the 
cloud resolving simulations. The summary is given in section 5. 

2. Model and experiments 

The cloud resolving model was originally developed by Soong and Ogura (1980), 
Soong and Tao (1980), and Tao and Simpson (1993). The 2D version of the model used 
by Sui et ai. (1994, 1998) and further modified by Li et al. (1999) is used in this study. 
The governing equations and model setup can be found in Li et al. (1999, 2002c). Several 
2D cloud resolving models have successfully simulated atmospheric thermodynamic 
profiles, cloud properties, and precipitation in the tropics during the Global Atmospheric 
Research Program Atlantic Tropical Experiment (GATE) (e.g., Xu and Randall, 1996; 
Grabowski et al., 1996; 1998; 1999) and Tropical Ocean Global Atmosphere Coupled 
Ocean- Atmosphere Response Experiment (TOGA CO ARE) (e.g., Wu et al., 1998; 1999; 
Li et al., 1999; 2002a;b;c; Johnson et al., 2004). It should be pointed out that the high- 
frequency temporal variability with spatial distributions such as individual cloud and 
associated dynamic and thermodynamic patterns might be distorted due to the 2D model 
setup. However, Grabowski et al. (1998) compared 3D and 2D simulations and showed 
the similar evolution in their thermodynamic fields, surface heat fluxes, and surface 
precipitation. 

For zonal mean heat budget as shown in Li et al. (1999), 
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( 1 ) 


dT Q m Q r it d(pw'e') _-o dd -o dT° 

— -j -rzr JX VU U . 

dt c p c p p dz dz dx 

Thus, the local change of zonal mean temperature is contributed by condensational 
heating, radiative heating, convergence of vertical heat flux, vertical temperature 
advection, and imposed horizontal temperature advection. 

For zonal mean moisture budget, 

( 2 ) 

dt p dz dz dx 

Thus, the local change of zonal mean water vapor is contributed by net condensation, 
convergence of vertical moisture flux, vertical moisture advection, and imposed 
horizontal moisture advection. 

The model is forced by the zonally uniform vertical velocity, zonal wind, and 
thermal and moisture advections, which are derived by Professor M. Zhang and his 
research group at the State University of New York at Stony Brook, based on the 6- 
hourly TOGA COARE observations within the Intensive Flux Array (IF A) region 
(Zhang, personal communication, 1999). The calculations are based on the constrained 
variational method on column-integrated budgets of mass, heat, moisture and momentum 
proposed by Zhang and Lin (1997). Hourly sea surface temperature at the Improved 
Meteorological (IMET) surface mooring buoy (1.75°S, 156°E) (Weller and Anderson, 
1996) is also imposed in the model. The model is integrated from 0400 LST 19 
December 1992 to 0400 LST 29 December 1992 (10 days total). Figure 1 shows the time 
evolution of vertical distribution of the large-scale atmospheric vertical velocity, zonal 
wind, and the time series of the sea surface temperature (SST) during the 10-day period. 
In this model setup, the horizontal boundary is periodic. The horizontal domain is 768 
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km, and the horizontal grid mesh is 1.5 km. The vertical grid resolution ranges from 
about 200 m near the surface to about 1 km about 100 mb. The time step is 12 s. 

The control experiment C uses hydrometeor equations (A1)-(A5) in the appendix, 
consisting of the full set of microphysical schemes used in the cloud resolving model. To 
test the impacts of the excluded microphysical processes in ZC97, we first identify the 14 
microphysical processes (the italicized terms in Fig. Al) being excluded by ZC97 in 
terms of microphysical schemes used in the cloud resolving model. A complete 
description of the terms can be found in the appendix (also see Li et al., 1999; 2002c). 
Then, the second experiment, C14, is carried out using (A10)-(A13) in which the 14 
microphysical terms are excluded. Note that (A10)-(A13) is derived for C14 that is 
symbolically similar to those used in ZC97, actually separate prognostic equations for the 
five hydrometeor species are still used in the simulation. 

We carry out the third experiment, CSFI, which is the same as C14 except for the 
inclusion of the P S fi(T<T 0 ). Thus, the CSFI uses (A10)-(A13) except that (A13b) is 
changed to 

P„ -PsAvr(.T<T.) + P„ a (T<r.) + P sr ,(T<T,). (3) 

The Psfi is the main process to consume cloud ice in the growth of snow. Hsie et 
al. (1980) modified the work of Orville and Kopp (1977) that was based on the equation 
of the rate of growth of ice crystals by deposition proposed by Koenig (1971), and 
formulated Psfi by the mixing ratio divided the timescale that is needed for an ice crystal 
to grow from radius 40pm to radius 50pm. Based on the aircraft observations, Krueger et 
al. (1995) suggested that the timescale in P S n should be for a crystal to grow from 40pm 
to radius 100pm, which increases mixing ratio of cloud ice as indicated in Li et al. 
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(1999). In this study, modified formulation of Psn by Krueger et al. (1995) is used. 
Krueger et al. (1995) called Psn as the snow formation associated with the Bergeron 
process (snow production from cloud ice via the growth of Bergeron-process embryos). 
However, Bergeron (Bergeron-Findeisen) process is the diffusional growth of ice crystals 
in the presence of supercooled water droplets. Therefore, the P S n is defined in this study 
as the depositional growth of snow from cloud ice. 

ZC97 calculated precipitation diagnostically based on the assumption that the 
precipitation water (q r ) and ice (q s and qg) fall out of clouds to ground immediately. A 
major effect being neglected by this as sumption is the interaction of precipitating 
particles with radiation. To examine the impacts of the precipitation-radiation interaction 
on tropical thermodynamics and cloud properties, the fourth experiment CN is carried 
out. The difference between the CN and C is that the CN excludes the precipitation- 
radiation interaction by simply setting zero mixing ratios of precipitation water and ice 
when the radiation is calculated. 

3. Important sink of cloud ice: depositional growth of snow 

We first examine the differences between experiments C and C14. Figure 2 shows 
the time series of surface rain rate and vertically integrated mass- weighted mixing ratios 
of cloud water ([qc]) and cloud ice ([qi])- The surface rain rates in the two experiments 
have similar temporal evolution because of the same imposed vertical velocity. However, 
[qc] and [qi] are much larger in C14 than in C, in particular, [qi] in C14 is more than one 
order of magnitude larger than in C (note the different plotting scales of [qi] for C 14 and 
C in Fig. 2c). 
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Since C14 has similar microphysical schemes as does ZC97, we wonder whether 
ZC97 has the same tendency to produce unrealistically large mixing ratio of cloud ice. 
However, the observational and GDAS data are not available during the TOGA COARE. 
Thus, we proceed to compare the ice water path (TWP: mass-integrated mixing ratio of 
ice clouds including cloud ice, snow and graupel) and liquid water path (LWP: mass- 
integrated mixing ratio of water clouds including cloud water and raindrops) simulated by 
the GDAS with those retrieved by NOAA/National Environmental Satellite, Data, 
Information Service (NESDIS)/Microwave Surface and Precipitation Products System 
(MSPPS) using the data from Advanced Microwave Sounder Unit (AJVISU) onboard 
NOAA-15, 16, and 17 in recent years. The LWP and IWP are retrieved using two 
AMSU-A window channels at 23.8 and 31.4 GHz (Weng et al., 1997; 2003) and two 
AMSU-B channels at 89 and 150 GHz (Weng and Grody, 2000; Weng et al., 2003), 
respectively. The horizontal resolution of the data used here is 50 km. 

Figure 3 shows the IWP versus LWP using the MSPPS and GDAS data over 
global oceanic tropics between 5°S and 5°N on March 2003 . The IWP simulated by the 
GDAS has the similar magnitudes to that observed in the MSPPS data (2 mm) whereas 
the LWP in the MSPPS data (3 mm) is three times as large as the GDAS simulated LWP 
(1 mm ) Since the IWP and LWP shown in Fig. 3a include both non-precipitating ice and 
water and precipitating ice and water whereas those shown in Fig. 3b include non- 
precipitation ice and water (cloud ice and cloud water) only, the comparison may suggest 
that the GDAS simulates the unrealistically large mixing ratio of ice clouds by using 
ZC97. Figures 2 and 3 also show that the cloud resolving model may simulate larger IWP 
and IWP in comparison with MSPPS observations. 
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The anomalous [qi] causes a large cooling above the ice clouds due to the 
reflection of the solar radiation and a large warming below due to the energy trapping of 
the longwave radiation, and positive difference of water vapor fields between C14 and C 
(Fig. 4). Table 1 also shows that 10-day mean fractional covers of raining convective, 
raining stratiform, non-raining stratiform clouds in C14 (C) are about 7.1%, 15.6%, and 
65.6%, (4.8%, 12.8%, and 52.4%) respectively. The result indicates a respective increase 
of 48%, 22%, and 25% in the fractional covers of raining convective, ra inin g stratiform, 
non-raining stratiform clouds from C to Cl 4. The C14 confirms unrealistic simulation of 
cloud ice. The budgets of [qi] in C and C14 are then analyzed to identify the cause of the 
unrealistic [qi] simulation in C14. 

The time series of vertically integrated budgets of [qi] (Eq. A3) in C and C14 is 
shown, respectively, in Fig. 5a and 5b. In C, the growth rate of [qi] by the vapor 
deposition ([Pdep]) is nearly balanced by the growth rate of snow by the conversion of 
cloud ice ([Psaut]) as well as by [P S fi] (Fig. 5a). Thus, the sum of the rates responsible 
for the growth of [qi] ([Sqi]) is very small. [Psfi] is smaller than [Psaut], but it is 
important sink of cloud ice. In the first five days of the integration of Cl 4, [Pdep] and 
[Psaut] cancel in large part, but there is a considerable magnitude of [Sqi] contributed by 
homogeneous freezing of cloud water ([Pihom]), which leads to the initial anomalous 
growth of cloud ice (Fig. 5b). The rapid decrease of the magnitude of [Psaut] cannot 
offset [Pdep] so that large [Sqi] is responsible for the rapid increase of [qi] in late evening 
of 24 December and early mo rnin g of 25 December (also see Fig. 3c). In the second five 
days of the integration, [Pihom] and [Psaci] (the collection of cloud ice by snow) become 
as important as [Pdep] and [Psaut] in the budget of [qi]. A comparison of budgets of [qi] 
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between C and C14 indicates that [Psfi] is important sink of cloud ice in C, but it is 
excluded in Cl 4, which may be the cause of anomalous growth of cloud ice in Cl 4. 

Thus, Psfi is included in the additional experiment CSFI that is conducted to further 
examine the cause of the difference between C14 and C. 

The temperature differences in the upper troposphere are sigmficantly reduced in 
CSFI — C (Fig. 6a) compared to those in C14— C (Fig. 4a). Also notice that the positive 
and negative temperature differences appear alternatively (Fig. 6a). The differences of 
water vapor ex hi bit a slightly positive trend with the most of the positive differences of 
less than 1 gkg' 1 (Fig. 6b). The fractional cloud covers in CSFI also display significant 
reduction from C14 and their 10-day mean values in CSFI become similar to those in C 
(Table 1). Thus, The CSFI confirms the crucial role of Psfi as the sink of [qi] in the 
balance of cloud ice. The results suggest that the budget of cloud ice be carefully treated 
to avoid unrealistic cloud calculations. 

Finally, 10-day mean vertically integrated cloud microphysics budgets in C and 
CSFI (Fig. 7) are compared. As shown in Li et al, (2002c), most of the conversion of 
cloud water to precipitation in C occurs primarily by the collection of cloud water by 
raindrops ([Pracw]=0.405 mmh' 1 ) at the temperature warmer than 0°C and by the riming 
of cloud water onto precipitation ice (snow and graupel) ([Pgacw]+ [Psacw]-0.157 mmh 
J ) at colder temperatures. Because the growth of precipitation ice through the conversion 
from cloud water is not allowed in CSFI, most of the conversion of cloud water to 
precipitation occurs p rimari ly by the collection of cloud water by raindrops 
([Pracw]=0.56 mmh ' 1 ) at the temperature warmer than 0°C. This leads to smaller mixing 
ratio of graupel in CSFI (0.07 mm) than in C (0.1 1 mm). In the budget of [qr], the 


evaporation of raindrops ([Prevp]=0.289 mmh' 1 ) is nearly compensated by the melting of 
graupel ([Pgmlt]-0.308 mmh' 1 ) so that the collection of cloud water by raindrops 
([PjEiACw]— 0.405 mmh ' 1 ) accounts mainly for surface rain rate (0.476 mmh' 1 ) in C. In 
CSFI, only half of [Prevp] (0.296 mmh' 1 ) is compensated by [Pgmlt] (0. 153 mmh" 1 ). In 
the budget of [qi], the vapor deposition rate in CSFI ([Pdep]=0.173 mmh 1 ) is larger than 
in C (0.124 mmh" 1 ). The rates of [Psaut], [Psfi], and [Psaci] in CSFI (0.126, 0.049, 0.006 
mmh" 1 ) are larger than in C (0.084, 0.035, 0.004 mmh' 1 ), respectively. The larger 
conversion rates from [qi] to [qs] in CSFI maintain the similar amount of [qs] of C. 

4. Important impacts on thermodynamics: precipitation-radiation interaction 

As mentioned already, ZC97 does not have the prognostic variables for 
precipitation (rain and precipitation ice), which leads to exclude the interaction of 
precipitating clouds with radiation in the model calculations in the GDAS. The effect of 
such simplification is examined by comparing CN and C. The comparison of 10-day 
mean fractional cloud covers between CN and C (Table 1) shows that the total fractional 
cloud covers in CN (71.5%) and C (70%) are slightly different. Then, a comparison study 
between CN and C is carried out to examine impacts of interaction between precipitation 
water and ice (raindrops, snow, and graupel) and radiation on atmospheric 
thermodynamics. 

Figure 8a shows the temporal and vertical distribution of temperature differences 
between CN and C (CN-C). Within the first three days of integration, the temperature 
differences are mostly positive except in late 21 December 1992 when negative 
temperature differences occur. From late 22 December 1992 on, the temperature 
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differences below 200 mb level become persistently negative with minimum values of 
—1.5 — 2°C. This indicates that the exclusion of precipitation-radiation interaction causes 
a cold bias in the troposphere. Li et al. (1999) carried out the comparison study between 
simulations and observations and showed that the simulated temperature in C is colder 
than the observed temperature in TOGA CO ARE (their Fig. 2a), the exclusion of the 
interaction processes in the simulation induces more cooling bias. 

To examine the physical processes responsible for the temperature differences, we 
integrated each term of (8) with time for CN and C separately and then took the 
difference (CN-C) to obtain the corresponding temperature differences due to 
condensational heating, radiative heating, convergence of vertical heat flux, and vertical 
temperature advection, respectively (Figs. 8b-8e). Note that the imposed horizontal 
temperature advection has no contribution to the temperature differences. The 
temperature differences due to radiation are negative below 500 mb, and positive above, 
with the maximum and minimum values reaching 6°C around 325 mb and — 6°C around 
575 mb, respectively, at the end of the integrations (Fig. 8c). The radiation-induced 
differences in thermal stratification cause a more stable layer above 500 mb in CN than in 
C. This corresponds to less ice hydrometeors in CN than in C (also see Fig. 9b), which 
lead to less latent heating above 500 mb in CN from 23 December 1992 on (Fig. 8b). The 
temperature differences due to the convergence of vertical heat flux in Fig. 8d show a 
banded structure with negative zones around 650-800 mb and 550 mb and mostly 
positive values elsewhere. The temperature differences due to vertical advection as 
shown in Fig. 8e are all negative except the 500-800 mb layer where positive values 
exist. Considering the contribution to the temperature difference of CN minus C by all 



terms, the negative temperature differences above 500 mb level are mainly contributed by 
less cloud heating as a result of more stability in CN, whereas more radiative cooling in 
CN is directly responsible for the negative temperature differences below 500 mb level. 

Next, we examine the temporal and vertical distribution of the difference fields of 
specific humidity (CN-C) in Fig. 9a, and the contribution to the moisture differences due 
to condensation, convergence of vertical moisture flux, and vertical moisture advection 
[see (9)] in Figs. 9b-d, respectively. The difference values of the water vapor mixing ratio 
are persistently negative below 500 mb during the integrations, with minimum around —1 
to —1.5 gkg' 1 . The significant drying begins on 23 December 1992. The exclusion of 
precipitation-radiation interaction in CN causes more unstable lower troposphere that 
leads to less evaporation of raindrops associated with subsidence that results in more 
drying, compared to C (Fig. 9b). Note that the imposed horizontal moisture advection has 
no contribution to the moisture differences. The contributions to vertical distributions of 
moisture differences by the other two terms, the convergence of vertical moisture flux 
and vertical moisture advection, largely cancel each other due to an out-of-phase relation 
(Figs. 9c, 9d). Thus, the less rain evaporation associated with large-scale subsidence as a 
result of more instability in the lower troposphere increases water vapor more slowly that 
leads to negative vapor differences. In summary, the comparison between CN and C 
indicates that exclusion of the interaction of precipitation with radiation in the 
simulations causes a significant cooling and drying bias in troposphere. 


5. Summary 
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We carry out a sensitivity study to examine the effect of certain neglected 
microphysical processes in the prognostic cloud scheme proposed by Zhao and Carr 
(1997) that has been used in operational numerical weather prediction models. The ZC97 
scheme excludes some microphysical processes and precipitation-radiation interaction for 
computational efficiency. Evaluating the effect of such simplifications in the large-scale 
models is not a simple task due to the interactions among different physical and dynamic 
processes. A lack of detailed microphysics measurements makes the task even more 
daunting. In this study, a cloud resolving model and relevant observations are used to 
examine the possible effects of the simplifications in the thermodynamic and cloud 
simulations. 

The sensitivity tests are carried out with a 2D cloud resolving model. The cloud 
model is integrated for 10 days with the imposed large-scale vertical velocity, zonal 
wind, and large-scale horizontal thermal and moisture advection observed and derived 
from the Tropical Ocean Global Atmosphere Coupled Ocean— Atmosphere Response 
Experiment (TOGA COARE). The experiment with the prognostic cloud scheme 
excluding the depositional growth of snow from cloud ice shows that the mixing ratio of 
cloud ice grows rapidly, which leads to more than 20% increase of fractional cloud 
covers and unrealistic vertical stratification compared to the experiment with full cloud 
microphysical parameterization package, indicating that the exclusion of the depositional 
growth of snow from cloud ice in the simulation could cause anomalous growth of cloud 
ice. Such a unrealistically large cloud ice simulation is also found in the GDAS data 
when they are compared with the satellite-retrieved cloud products. 
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The comparison between the experiments including and excluding precipitation- 
radiation interaction shows that the exclusion of the interaction in the model causes a 
coo lin g and drying bias. A further analysis of heat budgets reveals that the experiment 
excluding the interaction exhibits a more stable upper troposphere (above 500 mb) and a 
more unstable lower troposphere (below 500 mb) relative to the experiment including the 
interaction. The more stable upper troposphere suppresses the development of ice clouds 
that is responsible for the cooling bias whereas the more radiative cooling accounts 
directly for cooling bias in mid and lower troposphere in the experiment excluding the 
interaction. The analysis of moisture budgets shows that the suppression of rain 
evaporation as a result of less stable mid and lower troposphere induces the drying bias 
when the model excludes the interaction process. It should be notice that the schemes 
used in the cloud resolving model are different from those used in ZC97 even if the same 
microphysical processes are included Nevertheless, the study suggests that the 
depositional growth of snow from cloud ice and precipitating hydrometeor-radiation 
interaction be included in prognostic cloud schemes such as ZC97 for better prediction of 
atmospheric thermodynamics and cloud properties. 
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Table 1 10-day mean fractional covers (%) of raining convective (fcrcc), raining 
stratiform (fcrsc), non-raining stratiform (fcnrsc) clouds. 



fcrcc 

fcrsc 

fcnrsc 

c 

4.8 

12.8 

52.4 

C14 

7.1 

15.6 

65.6 

CSFI 

5.1 

13.3 

54.4 

CN 

4.5 

11.9 

55.1 
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Appendix 

Prognostic Equations of Hydrometeors and Microphysical Processes Parameterized 

in the Cloud Resolving Model 

The prognostic equations of hydrometeors in the cloud resolving model can be 
expressed as 
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where the microphysical processes in the terms of the right-hand side of (A1-A5) and 
corresponding schemes are listed in Table Al, and definitions and values of the other 
quantities can be found in Li et al. (1999, 2002c). 

To derive prognostic cloud scheme similar to Zhao and Carr (1997) in the frame 
of microphysics schemes in the cloud model, (Al) and (A3), and (A4) and (A5) are 
added, respectively, 

d ^ c t ? A = — —u(q c + q t ) -=-^-pw(q c + q t ) + P^ + P DEP -P r ~P sg ~ C{q c ,P sg ), 
dt dx p dz 


= - - u(q, * q „ ) - i - p O - w„ )(?, + q, ) - C(P r . P„ ) - C. {P r , P„ ) 


+ P„+C(,q„PJ + 0.-S,)P !I , er (T<TJ*a-S,)PoDEp<T<TJ-P m . TS (T>TJ 

~P^(T>T,\ (A7) 

and the (A2) can be rewritten as 

— - i— p(w - Wjy )q r + P r + C(P r , P sg ) + C m (P r ,P sg )~ P rrvp , (A8) 

dt dx p dz 


where 


P r ~ PsACW ( T>T o) + PrAVT + PrACW + PgACW (P > To )> 


?n 


P„-P i <u r (T<TJ + Pz, cr (T<T p ) + P m (T<TJ*P ua (T<T 0 ) + P 0Aa (T<T P ), (A9b) 
Ciyc’Psg) ~ PgACW^X < t>) ~ PsFW (P < P 0 )~ P GACW P < P'> 9 (A9c) 

C(P r , P„ ) - P „ 0 (2- > T„ ) - P ua (T < T p ) - P aAa , (J < T.) - P^ (T < T, ) - P CFR (T < T c ), 

(A9d) 

C, OP, , P„ ) - P^r CT > Tg ) + jw (r > r„). (A9e) 

P r and P sg are the sums of microphysics processes directly linking to the growth of 
precipitating water (q r ) and precipitating ice (q s and q*) respectively; C(q c ,P S g) is 
conversion term between cloud water and precipitating ice; C(P r ,P sg ) and Cm(P r ,P S g) are 
conversion terms between precipitating water and ice. 

In addition to the assumption that the precipitating hydrometeors fall immediately 
out of clouds without any interaction with radiation, Zhao and Carr’s (1997) prognostic 
cloud scheme excludes 14 terms (see italicized terms in Fig. Al; Also see Fig. 3 in ZC97 
including C(q c ,P sg ), C(P r ,P sg ), Psdep(T<T 0 ), Pgdep(T<T 0 ), Pgacw(T>T 0 ), Psfi(T<T 0 ), 
Pgaci(T<T 0 ), and Praci(T<T 0 ). Thus, (A6)-(A8) become 


a( V g,) - “(9, + 9, ) - ir- + «< ) + - P ' - p * • 

at ox p dz 

diQs+Qg) d Id— v \ r- rp p \ p 

— — = - — «(?,+ q s ) - =■— P( W ~ W TV X?, + ^g)~ C m ( P r > P sg ) ~ P sg 

dt dx p dz 

~ P MLTS (T>T a ) 

P MLTG (T>T a ), 


ggr 

dt 


d(uq r ) 

dx 


+ P r + C m ( P r > P sg )- P REVr > 

p dz 


(A10) 


(All) 
(A 12) 


where 

K-Pm^(T>T.)*P unr *P u cw, ( A1 3a ) 


?1 



P, t ~P s <m(T<T') + P SAa (T<T'), (A13b) 

CJP'.P^-P^AT > TJ + P gult (T > T.). (A13c) 

Note that the prognostic equations of precipitating hydrometeors (raindrops, snow, and 
graupel) are kept here since the cloud schemes associated with the growth of non- 
precipitating hydrometeors (cloud water and cloud ice) due to the conversions from 
precipitating hydrometeors (raindrops, snow, and graupel) are proportional to covariance 
between the mixin g ratios of precipitation and non-precipitating hydrometeors in cloud 
microphysics parameterization schemes used in this study. Also note that separate 
prognostic equations for five hydrometeor spices are used in this study though (A10)- 
(A 13) are symbolically derived 


9? 


Table A1 List of microphysical processes and their parameterization schemes in 
appendix. The schemes are Rutledge and Hobbs (1983, 1984; RH83, RH84), Lin et al. 
(1983, LFO), Tao et al. (1989, TSM), and Krueger et al. (1995, KFLC). 


Notation 

Description 

Scheme 

Pmltg 

Growth of vapor by evaporation of liquid from graupel surface 

RJH84 

Pmlts 

Growth of vapor by evaporation of melting snow 

RH83 

Prevp 

Growth of vapor by evaporation of raindrops 

RH83 

PlMLT 

Growth of cloud water by melting of cloud ice 

RH83 

P CND 

Growth of cloud water by the condensation of supersaturated 
vapor 

TSM 

Pgmlt 

Growth of raindrops by melting of graupel 

RH84 

PsMLT 

Growth of raindrops by melting of snow 

RH83 

Praci 

Growth of raindrops by the accretion of cloud ice 

RH84 

Pracw 

Growth of raindrops by the collection of cloud water 

RH83 

Pracs 

Growth of raindrops by the accretion of snow 

RH84 

Praut 

Growth of raindrops by the autoconversion of cloud water 

LFO 

PlDW 

Growth of cloud ice by the deposition of cloud water 

KFLC 

PlACR 

Growth of cloud ice by the accretion of rain 

RH84 

PlHOM 

Growth of cloud ice by the homogeneous freezing of cloud 
water 


Pdep 

Growth of cloud ice by the deposition of supersaturated vapor 

TSM 

PSAUT 

Growth of snow by the conversion of cloud ice 

RH83 

Psaci 

Growth of snow by the collection of cloud ice 

RH83 

PSACW 

Growth of snow by the accretion of cloud water 

RH83 

PsFW 

Growth of snow by the deposition and riming of cloud water 

KFLC 

PsFI 

Depositional growth of snow from cloud ice 

KFLC 

PsACR 

Growth of snow by the accretion of raindrops 

LFO 

PsDEP 

Growth of snow by the deposition of vapor 

RH83 

Pgaq 

Growth of graupel by the collection of cloud ice 

RH84 

Pgacr 

Growth of graupel by the accretion of raindrops 

RH84 

Pgacs 

Growth of graupel by the accretion of snow 

RH84 

Pgacw 

Growth of graupel by the accretion of cloud water 

RH84 

PwACS 

Growth of graupel by the riming of snow 

RH84 

Pgdep 

Growth of graupel by the deposition of vapor 

RH84 

Pgfr 

Growth of graupel by the freezing of raindrops 

LFO 
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Figure Captions 

Fig. 1 Temporal and vertical distribution of vertical velocity (a), zonal wind (b), and the 
time series of sea surface temperature (c) taken from TOGA COAJRE for 10-day 
period. Upward motion in (a) and westerly wind in (b) are shaded. Units of 
vertical velocity, zonal wind, and sea surface temperature are cm s' 1 , m s' 1 , and 
°C, respectively. 

Fig. 2 Time series of surface rain rate, Ps (a), vertically integrated mixing ratios of cloud 
water, [qc] (b), and cloud ice, [qi] (c) simulated in C (dashed) and Cl 4 (solid). 
The plotting scales for [qi] are 0-0.4 mm for C and 0-3 mm for C14 respectively. 
Units of Ps, [qc], [qi] are mm h' 1 , mm, mm. 

Fig. 3 IWP versus LWP using the MSPPS data (a) and IWP (cloud ice) versus LWP 
(cloud water) using the GDAS data (b). Unit is mm. 

Fig. 4 Temporal and vertical distribution of difference fields of (a) temperature (°C), and 
(b) specific humidity (gkg 1 ) for C14-C. Positive differences are shaded. 

Fig. 5 Time series of cloud ice budgets simulated in (a) C and (b) C14 respectively. All 
terms are vertically integrated. Unit is mm h' 1 . 

Fig. 6 Same as in Fig. 4 except for CSFI-C. 

Fig. 7 10-day mean vertically integrated cloud microphysics budgets simulated in (a) C 
and (b) CSFI. Units for cloud hydrometeors and conversions are mm and mm h" 1 , 
respectively. 

Fig. 8 Temporal and vertical distribution of (a) temperature differences between CN and 
C (CN-C), and differences due to (b) condensational heating, (c) radiative 
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heating, (d) convergence of vertical heat flux, and (e) vertical temperature 
advection. Unit is °C. Positive differences are shaded. 

Fig. 9 Temporal and vertical distribution of (a) differences of specific humidity between 
CN and C (CN-C), and the differences due to (b) condensation, (c) convergence 
of vertical moisture flux, and (d) vertical moisture advection. Unit is gkg' 1 . 
Positive differences are shaded. 

Fig. A1 Diagram of microphysical processes in the cloud resolving model. The terms in 
italics are excluded in the experiment C14. 







1993 


Fig. 1 Temporal and vertical distribution of vertical velocity (a), zonal wind (b), and the 
time series of sea surface temperature (c) taken from TOGA COARE for 10-day period. 
Upward motion in (a) and westerly wind in (b) are shaded. Units of vertical velocity, 
zonal wind, and sea surface temperature are cm s' 1 , m s' 1 , and °G, respectively. 
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Fig. 2 Time series of surface rain rate, Ps (a), vertically integrated mixing ratios of cloud 
water, [qc] (b), and cloud ice, [qi] (c) simulated in C (dashed) and Cl 4 (solid). The 
plotting scales for [qi] are 0-0.4 mm for C and 0-3 mm for C14 respectively. Units of Ps, 
[qc], [qi] are mm h' 1 , mm, mm. 
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Fig. 3 IWP versus LWP using the MSPPS data (a) and IWP (cloud ice) versus LWP 
(cloud water) using the GDAS data (b). Unit is mm. 
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Fig. 4 Temporal and vertical distribution of difference fields of (a) temperature (°C) and 
(b) specific humidity (gle g -1 ) for C14-C. Positive differences are shaded. 
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Fig. 8 Temporal and vertical distribution of (a) temperature differences between CN and 
C (CN-C), and differences due to (b) condensational beating, (c) radiative heating, (d) 
convergence of vertical heat flux, and (e) vertical temperature advection. Unit is °C. 
Positive differences are shaded. 





Fig. 9 Temporal and vertical distribution of (a) differences of specific humidity between 
CN and C (CN-C), and the differences due to (b) condensation, (c) convergence of 
vertical moisture flux, and (d) vertical moisture advection. Unit is gkg \ Positive 
differences are shaded. 
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Fig. A1 Diagram of microphysical processes in the cloud resolving model. The terms 
italics are excluded in the experiment C14. 
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Tropical Convective Responses to Microphysical and Radiative Processes: 

A Sensitivity Study With a 2D Cloud Resolving Model 

Xiaofan Li, C.-H. Sui, K.-M. Lau, W.-K. Tao 
Meteorology and Atmospheric Physics 

Popular Summary 

Prognostic cloud schemes are increasingly used in weather and climate models in order to 
better treat cloud-radiation processes. Simplifications are often made in such schemes for 
computational efficiency, like the scheme being used in the National Centers for Environment 
Prediction models that excludes some microphysical processes and precipitation-radiation 
interaction. In this study, sensitivity tests with a 2D cloud resolving model are carried out to 
examine effects of the excluded microphysical processes and precipitation-radiation 
interaction on tropical thermodynamics and cloud properties. The model is integrated for 1 0 
days with the imposed vertical velocity derived from die Tropical Ocean Global Atmosphere 
Coupled Ocean-Atmosphere Response Experiment. The experiment excluding the 
depositional growth of snow from cloud ice shows anomalous growth of cloud ice and more 
than 20% increase of fractional cloud cover, indicating that the lack of the depositional snow 
growth causes unrealistically large mixing ratio of cloud ice. The experiment excluding the 
precipitation-radiation interaction displays a significant cooling and drying bias. The analysis 
of heat and moisture budgets shows that the simulation without the interaction produces more 
stable upper troposphere and more unstable mid and lower troposphere than does the 
simulation with the interaction. Thus, the suppressed growth of ice clouds in upper 
troposphere and stronger radiative cooling in mid and lower troposphere are responsible for 
the cooling bias, and less evaporation of rain associated with the large-scale subsidence 
induces the drying in mid and lower troposphere. 



